CRISPR/Cas9-edited ROS1 + non-small cell lung cancer cell lines highlight differential drug sensitivity in 2D vs 3D cultures while reflecting established resistance profiles.

INTRODUCTION: The study of resistance-causing mutations in oncogene-driven tumors is fundamental to guide clinical decisions. Several point mutations affecting the ROS1 kinase domain have been identified in the clinical setting, but their impact requires further exploration, particularly in improved pre-clinical models. Given the scarcity of solid pre-clinical models to approach rare cancer subtypes like ROS1 + NSCLC, CRISPR/Cas9 technology allows the introduction of mutations in patient-derived cell lines for which resistant variants are difficult to obtain due to the low prevalence of cases within the clinical setting. METHODS: In the SLC34A2-ROS1 rearranged NSCLC cell line HCC78, we knocked-in through CRISPR/Cas9 technology three ROS1 drug resistance-causing mutations: G2032R, L2026M and S1986Y. Such variants are located in different functional regions of the ROS1 kinase domain, thus conferring TKI resistance through distinct mechanisms. We then performed pharmacological assays in 2D and 3D to assess the cellular response of the mutant lines to crizotinib, entrectinib, lorlatinib, repotrectinib and ceritinib. In addition, immunoblotting assays were performed in 2D-treated cell lines to determine ROS1 phosphorylation and MAP kinase pathway activity. The area over the curve (AOC) defined by the normalized growth rate (NGR_fit) dose-response curves was the variable used to quantify the cellular response towards TKIs. RESULTS: Spheroids derived from ROS1G2032R cells were significantly more resistant to repotrectinib (AOC fold change = - 7.33), lorlatinib (AOC fold change = - 6.17), ceritinib (AOC fold change = - 2.8) and entrectinib (AOC fold change = - 2.02) than wild type cells. The same cells cultured as a monolayer reflected the inefficacy of crizotinib (AOC fold change = - 2.35), entrectinib (AOC fold change = - 2.44) and ceritinib (AOC fold change = - 2.12) in targeting the ROS1 G2032R mutation. ROS1L2026M cells showed also remarkable resistance both in monolayer and spheroid culture compared to wild type cells, particularly against repotrectinib (spheroid AOC fold change = - 2.19) and entrectinib (spheroid AOC fold change = - 1.98). ROS1S1986Y cells were resistant only towards crizotinib in 2D (AOC fold change = - 1.86). Overall, spheroids showed an increased TKI sensitivity compared to 2D cultures, where the impact of each mutation that confers TKI resistance could be clearly distinguished. Western blotting assays qualitatively reflected the patterns of response towards TKI observed in 2D culture through the levels of phosphorylated-ROS1. However, we observed a dose-response increase of phosphorylated-Erk1/2, suggesting the involvement of the MAPK pathway in the mediation of apoptosis in HCC78 cells. CONCLUSION: In this study we knock-in for the first time in a ROS1 + patient-derived cell line, three different known resistance-causing mutations using CRISPR/Cas9 in the endogenous translocated ROS1 alleles. Pharmacological assays performed in 2D and 3D cell culture revealed that spheroids are more sensitive to TKIs than cells cultured as a monolayer. This direct comparison between two culture systems could be done thanks to the implementation of normalized growth rates (NGR) to uniformly quantify drug response between 2D and 3D cell culture. Overall, this study presents the added value of using spheroids and positions lorlatinib and repotrectinib as the most effective TKIs against the studied ROS1 resistance point mutations.

Identification of a DLG3 stop mutation in the MRX20 family.

Here, we identified the causal mutation in the MRX20 family, one of the larger X-linked pedigrees that have been described in which no gene had been identified up till now. In 1995, the putative disease gene had been mapped to the pericentromeric region on the X chromosome, but no follow-up studies were performed. Here, whole exome sequencing (WES) on two affected and one unaffected family member revealed the c.195del/p.(Thr66ProfsTer55) mutation in the DLG3 gene (NM_021120.4) that segregated with the affected individuals in the family. DLG3 mutations have been consequently associated with intellectual disability and are a plausible explanation for the clinical abnormalities observed in this family. In addition, we identified two other variants co-segregating with the phenotype: a stop gain mutation in SSX1 (c.358G>T/p.(Glu120Ter)) (NM_001278691.2) and a nonsynonymous SNV in USP27X (c.56 A>G/p.(Gln19Arg)) (NM_001145073.3). RNA sequencing revealed 14 differentially expressed genes (p value < 0.1) in 7 affected males compared to 4 unaffected males of the family, including four genes known to be associated with neurological disorders. Thus, in this paper we identified the c.195del/p.(Thr66ProfsTer55) mutation in the DLG3 gene (NM_021120.4) as likely responsible for the phenotype observed in the MRX20 family.

Pre-clinical modelling of ROS1+ non-small cell lung cancer.

Non-small cell lung cancer (NSCLC) is a heterogeneous group of diseases which accounts for 80% of newly diagnosed lung cancers. In the previous decade, a new molecular subset of NSCLC patients (around 2%) harboring rearrangements of the c-ros oncogene 1 was defined. ROS1+ NSCLC is typically diagnosed in young, nonsmoker individuals presenting an adenocarcinoma histology. Patients can benefit from tyrosine kinase inhibitors (TKIs) such as crizotinib and entrectinib, compounds initially approved to treat ALK-, MET- or NTRK- rearranged malignancies respectively. Given the low prevalence of ROS1-rearranged tumors, the use of TKIs was authorized based on pre-clinical evidence using limited experimental models, followed by basket clinical trials. After initiating targeted therapy, disease relapse is reported in approximately 50% of cases as a result of the appearance of resistance mechanisms. The restricted availability of TKIs active against resistance events critically reduces the overall survival. In this review we discuss the pre-clinical ROS1+ NSCLC models developed up to date, highlighting their strengths and limitations with respect to the unmet clinical needs. By combining gene-editing tools and novel cell culture approaches, newly developed pre-clinical models will enhance the development of next-generation tyrosine kinase inhibitors that overcome resistant tumor cell subpopulations.

varAmpliCNV: analyzing variance of amplicons to detect CNVs in targeted NGS data.

MOTIVATION: Computational identification of copy number variants (CNVs) in sequencing data is a challenging task. Existing CNV-detection methods account for various sources of variation and perform different normalization strategies. However, their applicability and predictions are restricted to specific enrichment protocols. Here, we introduce a novel tool named varAmpliCNV, specifically designed for CNV-detection in amplicon-based targeted resequencing data (Haloplex™ enrichment protocol) in the absence of matched controls. VarAmpliCNV utilizes principal component analysis (PCA) and/or metric dimensional scaling (MDS) to control variances of amplicon associated read counts enabling effective detection of CNV signals. RESULTS: Performance of VarAmpliCNV was compared against three existing methods (ConVaDING, ONCOCNV and DECoN) on data of 167 samples run with an aortic aneurysm gene panel (n = 30), including 9 positive control samples. Additionally, we validated the performance on a large deafness gene panel (n = 145) run on 138 samples, containing 4 positive controls. VarAmpliCNV achieved higher sensitivity (100%) and specificity (99.78%) in comparison to competing methods. In addition, unsupervised clustering of CNV segments and visualization plots of amplicons spanning these regions are included as a downstream strategy to filter out false positives. AVAILABILITY AND IMPLEMENTATION: The tool is freely available through galaxy toolshed and at: https://hub.docker.com/r/cmgantwerpen/varamplicnv. Supplementary Data File S1: https://tinyurl.com/2yzswyhh; Supplementary Data File S2: https://tinyurl.com/ycyf2fb4. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Non-Mendelian inheritance patterns and extreme deviation rates of CGG repeats in autism.

As expansions of CGG short tandem repeats (STRs) are established as the genetic etiology of many neurodevelopmental disorders, we aimed to elucidate the inheritance patterns and role of CGG STRs in autism-spectrum disorder (ASD). By genotyping 6063 CGG STR loci in a large cohort of trios and quads with an ASD-affected proband, we determined an unprecedented rate of CGG repeat length deviation across a single generation. Although the concept of repeat length being linked to deviation rate was solidified, we show how shorter STRs display greater degrees of size variation. We observed that CGG STRs did not segregate by Mendelian principles but with a bias against longer repeats, which appeared to magnify as repeat length increased. Through logistic regression, we identified 19 genes that displayed significantly higher rates and degrees of CGG STR expansion within the ASD-affected probands (P < 1 × 10-5). This study not only highlights novel repeat expansions that may play a role in ASD but also reinforces the hypothesis that CGG STRs are specifically linked to human cognition.

Protein interaction network analysis reveals genetic enrichment of immune system genes in frontotemporal dementia.

Interactors of protein products of known genes for frontotemporal dementia (FTD) are likely to be involved in the molecular pathways towards disease. We therefore applied protein interaction network (PIN) analysis to prioritize candidate genes for rare variant association analysis. We created an FTD-PIN starting from known FTD genes downloading their physical interactors and performed functional enrichment analyses. We identified overrepresented processes in FTD and selected genes (n = 440) belonging to these processes for rare variant analysis in a Belgian cohort of 228 FTD patients and 345 controls. SKAT-O analysis suggested TNFAIP3 as the top gene (p = 0.7 × 10-3) reaching near test-wide significance (p = 2.5 × 10-4). We then analyzed the TNFAIP3-subnetwork within the FTD-PIN which indicated enrichment of several immune signaling networks, suggesting that disrupted immune signaling may be implicated in TNFAIP3-related FTD. Our study demonstrates that integration of PINs with genetic data is a useful approach to increase the power for rare variant association analysis. Furthermore, we present a computational pipeline for identifying potential novel therapeutic targets and risk-modifying variants.

Molecular characterization and investigation of the role of genetic variation in phenotypic variability and response to treatment in a large pediatric Marfan syndrome cohort.

PURPOSE: In a large cohort of 373 pediatric patients with Marfan syndrome (MFS) with a severe cardiovascular phenotype, we explored the proportion of patients with MFS with a pathogenic FBN1 variant and analyzed whether the type/location of FBN1 variants was associated with specific clinical characteristics and response to treatment. Patients were recruited on the basis of the following criteria: aortic root z-score > 3, age 6 months to 25 years, no prior or planned surgery, and aortic root diameter < 5 cm. METHODS: Targeted resequencing and deletion/duplication testing of FBN1 and related genes were performed. RESULTS: We identified (likely) pathogenic FBN1 variants in 91% of patients. Ectopia lentis was more frequent in patients with dominant-negative (DN) variants (61%) than in those with haploinsufficient variants (27%). For DN FBN1 variants, the prevalence of ectopia lentis was highest in the N-terminal region (84%) and lowest in the C-terminal region (17%). The association with a more severe cardiovascular phenotype was not restricted to DN variants in the neonatal FBN1 region (exon 25-33) but was also seen in the variants in exons 26 to 49. No difference in the therapeutic response was detected between genotypes. CONCLUSION: Important novel genotype-phenotype associations involving both cardiovascular and extra-cardiovascular manifestations were identified, and existing ones were confirmed. These findings have implications for prognostic counseling of families with MFS.

A Panel-Based Sequencing Analysis of Patients with Paget's Disease of Bone Suggests Enrichment of Rare Genetic Variation in regulators of NF-κB Signaling and Supports the Importance of the 7q33 Locus.

Paget's disease of bone (PDB) is a common bone disorder characterized by focal lesions caused by increased bone turnover. Monogenic forms of PDB and PDB-related phenotypes as well as genome-wide association studies strongly support the involvement of genetic variation in components of the NF-κB signaling pathway in the pathogenesis of PDB. In this study, we performed a panel-based mutation screening of 52 genes. Single variant association testing and a series of gene-based association tests were performed. The former revealed a novel association with NFKBIA and further supports an involvement of variation in NR4A1, VCP, TNFRSF11A, and NUP205. The latter indicated a trend for enrichment of rare genetic variation in GAB2 and PRKCI. Both single variant tests and gene-based tests highlighted two genes, NR4A1 and NUP205. In conclusion, our findings support the involvement of genetic variation in modulators of NF-κB signaling in PDB and confirm the association of previously associated genes with the pathogenesis of PDB.

Resequencing of candidate genes for Keratoconus reveals a role for Ehlers-Danlos Syndrome genes.

The involvement of genetic factors in the pathogenesis of KC has long been recognized but the identification of variants affecting the underlying protein functions has been challenging. In this study, we selected 34 candidate genes for KC based on previous whole-exome sequencing (WES) and the literature, and resequenced them in 745 KC patients and 810 ethnically matched controls from Belgium, France and Italy. Data analysis was performed using the single variant association test as well as gene-based mutation burden and variance components tests. In our study, we detected enrichment of genetic variation across multiple gene-based tests for the genes COL2A1, COL5A1, TNXB, and ZNF469. The top hit in the single variant association test was obtained for a common variant in the COL12A1 gene. These associations were consistently found across independent subpopulations. Interestingly, COL5A1, TNXB, ZNF469 and COL12A1 are all known Ehlers-Danlos Syndrome (EDS) genes. Though the co-occurrence of KC and EDS has been reported previously, this study is the first to demonstrate a consistent role of genetic variants in EDS genes in the etiology of KC. In conclusion, our data show a shared genetic etiology between KC and EDS, and clearly confirm the currently disputed role of ZNF469 in disease susceptibility for KC.

Abundancy of polymorphic CGG repeats in the human genome suggest a broad involvement in neurological disease.

Expanded CGG-repeats have been linked to neurodevelopmental and neurodegenerative disorders, including the fragile X syndrome and fragile X-associated tremor/ataxia syndrome (FXTAS). We hypothesized that as of yet uncharacterised CGG-repeat expansions within the genome contribute to human disease. To catalogue the CGG-repeats, 544 human whole genomes were analyzed. In total, 6101 unique CGG-repeats were detected of which more than 93% were highly variable in repeat length. Repeats with a median size of 12 repeat units or more were always polymorphic but shorter repeats were often polymorphic, suggesting a potential intergenerational instability of the CGG region even for repeats units with a median length of four or less. 410 of the CGG repeats were associated with known neurodevelopmental disease genes or with strong candidate genes. Based on their frequency and genomic location, CGG repeats may thus be a currently overlooked cause of human disease.

Delineation of a new fibrillin-2-opathy with evidence for a role of FBN2 in the pathogenesis of carpal tunnel syndrome.

BACKGROUND: Although carpal tunnel syndrome (CTS) is the most common form of peripheral entrapment neuropathy, its pathogenesis remains largely unknown. An estimated heritability index of 0.46 and an increased familial occurrence indicate that genetic factors must play a role in the pathogenesis. METHODS AND RESULTS: We report on a family in which CTS occurred in subsequent generations at an unusually young age. Additional clinical features included brachydactyly and short Achilles tendons resulting in toe walking in childhood. Using exome sequencing, we identified a heterozygous variant (c.5009T>G; p.Phe1670Cys) in the fibrillin-2 (FBN2) gene that co-segregated with the phenotype in the family. Functional assays showed that the missense variant impaired integrin-mediated cell adhesion and migration. Moreover, we observed an increased transforming growth factor-ß signalling and fibrosis in the carpal tissues of affected individuals. A variant burden test in a large cohort of patients with CTS revealed a significantly increased frequency of rare (6.7% vs 2.5%-3.4%, p<0.001) and high-impact (6.9% vs 2.7%, p<0.001) FBN2 variants in patient alleles compared with controls. CONCLUSION: The identification of a novel FBN2 variant (p.Phe1670Cys) in a unique family with early onset CTS, together with the observed increased frequency of rare and high-impact FBN2 variants in patients with sporadic CTS, strongly suggest a role of FBN2 in the pathogenesis of CTS.

Author Correction: Large-scale targeted sequencing identifies risk genes for neurodevelopmental disorders.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Large-scale targeted sequencing identifies risk genes for neurodevelopmental disorders.

Most genes associated with neurodevelopmental disorders (NDDs) were identified with an excess of de novo mutations (DNMs) but the significance in case-control mutation burden analysis is unestablished. Here, we sequence 63 genes in 16,294 NDD cases and an additional 62 genes in 6,211 NDD cases. By combining these with published data, we assess a total of 125 genes in over 16,000 NDD cases and compare the mutation burden to nonpsychiatric controls from ExAC. We identify 48 genes (25 newly reported) showing significant burden of ultra-rare (MAF < 0.01%) gene-disruptive mutations (FDR 5%), six of which reach family-wise error rate (FWER) significance (p < 1.25E-06). Among these 125 targeted genes, we also reevaluate DNM excess in 17,426 NDD trios with 6,499 new autism trios. We identify 90 genes enriched for DNMs (FDR 5%; e.g., GABRG2 and UIMC1); of which, 61 reach FWER significance (p < 3.64E-07; e.g., CASZ1). In addition to doubling the number of patients for many NDD risk genes, we present phenotype-genotype correlations for seven risk genes (CTCF, HNRNPU, KCNQ3, ZBTB18, TCF12, SPEN, and LEO1) based on this large-scale targeted sequencing effort.

Overrepresentation of genetic variation in the AnkyrinG interactome is related to a range of neurodevelopmental disorders.

Upon the discovery of numerous genes involved in the pathogenesis of neurodevelopmental disorders, several studies showed that a significant proportion of these genes converge on common pathways and protein networks. Here, we used a reversed approach, by screening the AnkyrinG protein-protein interaction network for genetic variation in a large cohort of 1009 cases with neurodevelopmental disorders. We identified a significant enrichment of de novo potentially disease-causing variants in this network, confirming that this protein network plays an important role in the emergence of several neurodevelopmental disorders.

Osmotic stress inhibits leaf growth of Arabidopsis thaliana by enhancing ARF-mediated auxin responses.

We investigated the interaction between osmotic stress and auxin signaling in leaf growth regulation. Therefore, we grew Arabidopsis thaliana seedlings on agar media supplemented with mannitol to impose osmotic stress and 1-naphthaleneacetic acid (NAA), a synthetic auxin. We performed kinematic analysis and flow-cytometry to quantify the effects on cell division and expansion in the first leaf pair, determined the effects on auxin homeostasis and response (DR5::ß-glucuronidase), performed a next-generation sequencing transcriptome analysis and investigated the response of auxin-related mutants. Mannitol inhibited cell division and expansion. NAA increased the effect of mannitol on cell division, but ameliorated its effect on expansion. In proliferating cells, NAA and mannitol increased free IAA concentrations at the cost of conjugated IAA and stimulated DR5 promotor activity. Transcriptome analysis shows a large overlap between NAA and osmotic stress-induced changes, including upregulation of auxin synthesis, conjugation, transport and TRANSPORT INHIBITOR RESPONSE1 (TIR1) and AUXIN RESPONSE FACTOR (ARF) response genes, but downregulation of Aux/IAA response inhibitors. Consistently, arf7/19 double mutant lack the growth response to auxin and show a significantly reduced sensitivity to osmotic stress. Our results show that osmotic stress inhibits cell division during leaf growth of A. thaliana at least partly by inducing the auxin transcriptional response.

Insufficient evidence for a role of SERPINF1 in otosclerosis.

Otosclerosis is a common form of hearing loss (HL) due to abnormal remodeling of the otic capsule. The genetic causes of otosclerosis remain largely unidentified. Only mutations in a single gene, SERPINF1, were previously published in patients with familial otosclerosis. To unravel the contribution of genetic variation in this gene to otosclerosis, this gene was re-sequenced in a large population of otosclerosis patients and controls. Resequencing of the 5' and 3' UTRs, coding regions, and exon-intron boundaries of SERPINF1 was performed in 1604 unrelated otosclerosis patients and 1538 unscreened controls, and in 62 large otosclerosis families. Our study showed no enrichment of rare variants, stratified by type, in SERPINF1 in patients versus controls. Furthermore, the c.392C > A (p.Ala131Asp) variant, previously reported as pathogenic, was identified in three patients and four controls, not replicating its pathogenic nature. We could also not find evidence for a pathogenic role in otosclerosis for 5' UTR variants in the SERPINF1-012 transcript (ENST00000573763), described as the major transcript in human stapes. Furthermore, no rare variants were identified in the otosclerosis families. This study does not support a pathogenic role for variants in SERPINF1 as a cause of otosclerosis. Therefore, the etiology of the disease remains largely unknown and will undoubtedly be the focus of future studies.

Copy number variation analysis in bicuspid aortic valve-related aortopathy identifies TBX20 as a contributing gene.

Bicuspid aortic valve (BAV) is the most common congenital heart defect (CHD), affecting 1-2% of the population. BAV is associated with thoracic aortic aneurysms (TAAs). Deleterious copy number variations (CNVs) were found previously in up to 10% of CHD cases. This study aimed at unravelling the contribution of deleterious deletions or duplications in 95 unrelated BAV/TAA patients. Seven unique or rare CNVs were validated, harbouring protein-coding genes with a role in the cardiovascular system. Based on the presence of overlapping CNVs in patients with cardiovascular phenotypes in the DECIPHER database, the identification of similar CNVs in whole-exome sequencing data of 67 BAV/TAA patients and suggested topological domain involvement from Hi-C data, supportive evidence was obtained for two genes (DGCR6 and TBX20) of the seven initially validated CNVs. A rare variant burden analysis using next-generation sequencing data from 637 BAV/TAA patients was performed for these two candidate genes. This revealed a suggestive genetic role for TBX20 in BAV/TAA aetiology, further reinforced by segregation of a rare TBX20 variant with the phenotype within a BAV/TAA family. To conclude, our results do not confirm a significant contribution for deleterious CNVs in BAV/TAA as only one potentially pathogenic CNV (1.05%) was identified. We cannot exclude the possibility that BAV/TAA is occasionally attributed to causal CNVs though, or that certain CNVs act as genetic risk factors by creating a sensitised background for BAV/TAA. Finally, accumulative evidence for TBX20 involvement in BAV/TAA aetiology underlines the importance of this transcription factor in cardiovascular disease.

Confirmation of the role of pathogenic SMAD6 variants in bicuspid aortic valve-related aortopathy.

Progressive dilatation of the thoracic aorta leads to thoracic aortic aneurysm (TAA), which is often asymptomatic but predisposes to lethal aortic dissections and ruptures. TAA is a common complication in patients with bicuspid aortic valve (BAV). Recently, rare loss-of-function SMAD6 variants were shown to contribute significantly to the genetic aetiology of BAV/TAA. Intriguingly, patients with craniosynostosis have also been reported to be explained molecularly by similar loss-of-function SMAD6 variants. While significantly reduced penetrance of craniosynostosis has been reported for the SMAD6 variants as such, near-complete penetrance is reached upon co-occurrence with a common BMP2 SNP risk allele. Here, we report on the results of a SMAD6-variant analysis in 473 unrelated non-syndromic TAA patients, of which the SMAD6-positive individuals were also studied for the presence of the BMP2 risk allele. Although only 14% of the TAA patients also presented BAV, all novel likely pathogenic SMAD6 variants (N = 7) were identified in BAV/TAA individuals, further establishing the role of SMAD6 variants to the aetiology of BAV/TAA and revealing limited contribution to TAA development in patients with a tricuspid aortic valve. Familial segregation studies confirmed reduced penetrance (82%) and variable clinical expressivity, with coarctation of the aorta being a common comorbidity. None of our six BMP2+/SMAD6+ patients presented with craniosynostosis. Hence, the proposed digenic model for craniosynostosis was not supported in the presented BAV/TAA cohort, suggesting that additional factors are at play. Finally, our data provide improved insights into the clinical spectrum of SMAD6-related BAV/TAA and has important implications for molecular diagnostics.

Genetic Variation in RIN3 in the Belgian Population Supports Its Involvement in the Pathogenesis of Paget's Disease of Bone and Modifies the Age of Onset.

Paget's disease of bone (PDB) is a common, late-onset bone disorder characterized by focal increase of bone turnover. Mutations in the SQSTM1 gene are found in up to 40% of patients and recent GWAS have led to novel associations with several loci. RIN3, the candidate gene located at the associated 14q32 locus, has recently been studied in a British cohort to elucidate its contribution to the pathogenesis. In this study, we performed a genetic screening of RIN3 in an unrelated cohort to validate these findings and to further explore genetic variation in this gene in the context of PDB. In our screening, we examined the 5' untranslated region (UTR), the exonic regions and the intron-exon boundaries of the gene in a control cohort and a patient cohort. Our findings show clustering of variation similar to the British cohort and support a protective role for common genetic variation (rs117068593, p.R279C) in the proline-rich region and a functionally relevant role for rare genetic variation in the domains that mediate binding and activation of its interaction partner, Rab5. Additive regression models, fitted for the common variants, validated the association of the rs117068593 variant with the disease (OR+/+ 0.315; OR+/- 0.562). In addition, our analyses revealed a potentially modifying effect of this variant on the age of onset of the disease. In conclusion, our findings support the involvement of genetic variation in RIN3 in PDB and suggest a role for RIN3 as a potential modifier of the age of onset of the disease.